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ABSTRACT. We study here the existence of solitary wave solutions of a generalized two- 
component Camassa-Holm system. In addition to those smooth solitary-wave solutions, 
we show that there are solitary waves with singularities: peaked and cusped solitary waves. 
We also demonstrate that all smooth solitary waves are orbitally stable in the energy space. 
We finally give a sufficient condition for global strong solutions to the equation in some 
special case. 

1. INTRODUCTION 

There are several classical models describing the motion of waves at the free surface of 
shallow water under the influence of gravity. Among these models, the best known is the 
Korteweg-de Vries (KdV) equation BT1I451 

Ut + Guu x + u xxx = 0. 

The KdV equation admits solitary wave solutions, i.e. solutions of the form u(t, x) — 
cp(x — ct) which travel with fixed speed c, and that vanish at infinity. The KdV solitary 
waves are smooth and retain their individuality under interaction and eventually emerge 
with their original shapes and speeds [20|. Moreover, KdV is an integrable infinite- 
dimensional Hamiltonian system (39). However, the KdV equation does not model the 
phenomenon of breaking for water waves. Instead, as soon as the initial profile uq e 
H 1 (M.), the solutions are global in time 11321 . whereas some shallow water waves break 

E3. 

Another model, the Camassa-Holm (CH) equation |3] 

Ut - u xxt + 3uu x — 2u x u xx + uu xxx (1.1) 

arises as a model for the unidirectional propagation of shallow water waves over a flat 
bottom ||3] [T2l [2T| |22, 23 30], as well as water waves moving over an underlying shear 
flow ||3T1 . Equation ((TTTJ is completely integrable with the Lax pair O and with infinitely 
many conservation laws as an bi-Hamiltonian system |25 1. 

The CH equation has many remarkable properties that KdV does not have like solitary 
waves with singularities and breaking waves. The CH equation admits peaked solitary 
waves or "peakons" |Q~| G] El : u(t,x) — ce~' x ~ ct ', c ^ 0, which are smooth except at 
the crests, where they are continuous, but have a jump discontinuity in the first derivative. 
The peakons capture a feature that is characteristic for the waves of great height - waves 
of the largest amplitude that are exact solutions of the governing equations for water waves 
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ll8l [T3ll44 ll. The CH equation also models wave breaking (i.e. the solution remains bounded 
while its slope becomes unbounded in finite time) ll3l 151 151171 [TOl [39l l45l . 

The CH equation also admits many multi -component generalizations. It is intriguing to 
know if the above two properties may persist in the systems. In this paper we consider the 
following generalized two-component CH system established in [4| which can be derived 
from shallow water theory with nonzero constant vorticity 

Ut - utxx - Au x + 3uu x - a(2u x u xx + uu xxx ) + pp x = 0, 
pt + (f>u)x = 0, 
or equivalently, using the linear momentum m = u — u xx , 

m t + <yum x — Au x + 2amu x + 3(1 — <j)uu x + pp x = 0, 
Pt + (pu) x = 0, 

where u(t, x) is the horizontal velocity and p(t, x) is related to the free surface elevation 
from equilibrium (or scalar density) with the boundary assumptions u — s- 0, p — > 1 as 
|x| — > oo. The scalar A > characterizes a linear underlying shear flow and hence system 
( |1.2| i models wave-current interactions. The real dimensionless constant a is a parameter 
which provides the competition, or balance, in fluid convection between nonlinear steep- 
ening and amplification due to stretching. 

When a = 1 it recovers the standard two-component CH system which is completely 
integrable [11 , 29, 42] as it can be written as compatibility condition of two linear systems 
(Lax pair) with a spectral parameter (, that is, 

-CV + c(m-^ + ^) *, 
In the case p = 0, it becomes 

+ 3uu x = a (2u x u xx + uu xxx ) , (1.3) 

which models finite length, small amplitude radial deformation waves in cylindrical hyper- 
elastic rods |19|. System \\.2\ has the following two Hamiltonians 

Hi = \ ( (u 2 + ul + (p - l) 2 ) dx, 

H 2 = - I (u 3 + auu\ + 2u(p - 1) + u(p - l) 2 - Au 2 ) dx. 

We study solitary wave solutions of \\.2\ , i.e. solutions of the form 

(u(x, t), p(x, t)) = (<p(x — ct), p(x — ct)), eel 

for some ip, p : K — > K such that <p — > 0, p — > 1 as \x\ — > oo. In the study of the CH 
traveling waves it was observed through phase-plane analysis B71 that both peaked and 
cusped traveling waves exist. Subsequently, Lenells l33l l34ll used a suitable framework 
for weak solutions to classify all weak traveling waves of the CH equation (JTTTJ and the 
hyperelastic rod equation ( |1.3| l. 

Using a natural weak formulation of the two-component CH system, we will establish 
exactly in what sense the peaked and cusped solitary waves are solutions. It was shown 
in [11, 38, 40 1 that when a = 1 the two component system \\.2\ has only smooth soli- 
tary waves, with a single crest profile and exponential decay far out. In [28 1, the authors 
considered a modified two-component CH equation which allows dependence on average 
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density as well as pointwise density and a linear dispersion is added to the first equation 
of the system. They showed that the modified system admits peakon solutions in both u 
and p. However it is unclear whether the generalized two-component CH system jl.2\ has 
solitary waves with singularities. We show here peaked solitary waves exist when a > 1. 
We also provide an implicit formula for the peaked solitary waves. However whether these 
peaked solitary waves are solitons still remains open. 

The stability of solitary waves is one of the fundamental qualitative properties of the 
solutions of nonlinear wave equations II43I . Due to the fact that the solitons hardly interact 
with each other at all it is reasonable to expect that they are stable. It has been proved 
that for the CH equation, the smooth solitary waves are orbitally stable ifTTIl . Moreover, the 
peakons, whether solitary waves or periodic waves, are also orbitally stable lfl5l[T6l[35l[36l . 
It was shown in [38 1 that when a = 1 all solitary waves are orbitally stable. We prove in 
this paper that when er < 1 all smooth solitary waves are orbitally stable. The proof of the 
stability basically follows the general approach in |26l . In comparison with the spectral 
arguments on the Hessian operator in |26l . here we require more precise analysis on the 
spectrum of a linearized operator around the solitary waves for the system \\.2\ . 

A special case of system ( |1.2| i is when a = 0. In the scalar equation case when a = it 
is the BBM equation |2|. The solutions are shown to be global in time. We show that the 
same results hold in the system case. 

This paper is organized as follows. In Section|2]we classify the solitary waves of ( [1.2) 1. 
In particular we show the existence of peaked solitary waves for a > 1. In Section [3] we 
prove that when a < 1 all smooth solitary waves are nonlinearly stable. Finally in Section 
|4] we show that the system ( |1.2[ ) is globally well-posed for a = 0. 



2. Solitary waves 

Let X — H 1 (M.) x L 2 (R) be a real Hilbert space with inner product (, ), and denote 
its element by u = {u,rf). The dual of X is X* = H^ 1 (R) x L 2 (R) and a natural 
isomorphism / from X to X* can be defined by 

l-d 2 x 
1 



/ = 



Using the map /, the paring (, ) between X and X* can be represented as 

{Iu,v) = (u,v) 1 + {r),£) , for u = (u, rj) € X, v= e X* , 

where (, ) s denotes the H S (R) x H~ S (W.) dual pairing. We will identify the second dual 
X** with X in a natural way. 

Since p — s- 1 as \x\ — > oo in ( |1.2) , we can let p = 1 + 77 with 77 — >• as \x\ — > 00 and 
hence we can rewrite system ( |1.2[ ) as 

u t - Utxx - Au x + 3uu x - a(2u x u xx + uu xxx ) + (1 + r})r\ x = 0, 
m + ((l + rj)u) x = 0. 

The two Hamiltonians introduced in the Introduction define the following two functionals 
onl 

E(u) = I I {u 2 +u 2 x + rj 2 ) dx, (2.2) 

F(u) = ~ I (u 3 + auu 2 x + 2w] + wq 2 - Au 2 ) dx, (2.3) 
2 7r 



with u = (n, 77) G X. The quantity E{u) associates with the translation invariance of 
\2.\\ . Using functional F(u), system ( |2. 1 ] > can be written in an abstract Hamiltonian form 

u t = JF'(u), (2.4) 

where J is a closed skew symmetric operator given by 



J = 



-d x (l-d 2 x )-l 






-dx 



and F' (u) : X ^> X* is the variational derivative of F in X at u. 
Note that if 

p(x) 



-e - lxl ,xe 



(2.5) 



then (1 — c^) 1 / = p * / for all / G L 2 (R). We can then further rewrite system ( |2.1[ ) in 
a weak form as 



"/ + "mi, + 'I V ' { -Au + 3 -y^" 2 + + ~(1 + j = 



(2.6) 



. % + ((1 + V)u) x = 0. 



Definition 2.1. L<?f < T < oo. A function u — {u,rf) G C([0,T);X) is called a 
solution of ( |2.1| l on [0, T) if if satisfies \2.6) in the distribution sense on [0, T) anj -E-(w) 
ana? f (it) are conserved. 

Applying transport equation theory combined with the method of Besov spaces, one 
may follow the similar argument as in ll27l to obtain the following local well-posedness 
result for the system ( |2.1| >. 

Theorem 2.2. // (uoi^o) G F[ s x _ff s_1 , s > 3/2, fnen fnere ewsf a maximal time 
T = ^(11 ( u o> 7 7o)||_fr s xH s - 1 ) > anci a unique solution (u, if) of ( |1.2| > in C([0, T); H s x 
i/ s_1 ) n ^([OjT);^* -1 x H s ~ 2 ) with (u,r))\ t = = («o,»7o)- Moreover, the solution 
depends continuously on the initial data and T is independent of s. 

It is easily seen from the embedding L°°(R) that -E(u) and are both 

well defined in H s x H s ~ l with s > 3/2, and E(u) is conserved, as suggested in the local 
wellposedness Thoerem |2.2| From \2A\ we see that 

d 



dt 



F{u) = (F'(u),u t ) = (F'(u),JF'(u)) = 0. 



So F(u) is also invariant. 

Now we give the definition of solitary waves of \2.l) . 

Definition 2.3. A solitary wave of ( |2.1| i is a nontrivial traveling wave solution of ( |2.1| i of 

the form (p c {t, x) = (ip c (x — ct), r\ c (x — ct)) G H 1 (M.) x with c € R anof </J c , ry c 

vanishing at infinity. 



For a solitary wave tp = ((p, 77) with speed c G R, it satisfies 



-cip 



■ p * 



. 3 ~ (To <7 9 1 „ 

-Ap + -g-V 2 + 2^ + 2 (1 + r/) 



-07 + (1 + 77)^ = 0, 



inP'(R). 
(2.7) 



Integrating the above system and applying (1 — d x ) to the first equation we get 

| -(c + A) V + cip xx + ^ 2 =a(ptp xx +°(p 2 x -^{l + T ] ) 2 + l i , inD'(M). 
{ —ci] + (1 + r))tp = 0. 

(2.8) 

The fact that the second equation of the above holds in a strong sense comes from the 
regularity of <p and 77. 

Proposition 2.4. // (ip, 77) is a solitary wave of ( |2.1[ ) for some c € R, then c ^ cmd 
^ cfor any i£i 

Proof. From the definition of solitary waves and the embedding theorem we know that ip 
and rj are both continuous. If c = 0, then (|2.8|> becomes 



{. 3 2 O" 2 1 . . o 1 

-Atp+~<p = crcpip xx + -<p x - -(1 + 77) +-, (29) 
(1 + ^ = 0. 

Since rj vanishes at infinity, the second equation of the above system indicates that tp(x) — 
for \x\ large enough. Denote xq — max{i : (f(x) ^ 0}. Hence <p(x) = on [xq, oo) 
and <p ^ on (xq — S, xq) for any 6 > 0. Consider now the first equation of ( |2.9| l on 
[xo, oo) we see that 77 = on [xq, 00). Then the continuity of 77 implies that there exists a 
<5i > such that 1 + rj(x) > on (xq — Si, Xq). This together with the second equation of 
( |2.9| > leads to <p(x) = on (xq — Si, Xq), which is a contradiction. Therefore c 7^ 0. 

Next we show (p ^ c. If not and there is some a?i € K such that <f{xi) = c. Then the 
second equation of ( |2.8| > infers that 

= (c - (f(xi)) j](xi) = 0, 
so c = 0, which is a contraction. □ 
Using the above proposition we obtain from the second equation of ( |2.8| l that 

»7=-^-. (2.10) 

c — ip 

Plugging this into the first equation of ( |2.8| l we obtain an equation for the unknown <p only 

- (c + A)<p + cp xx + ^ip 2 =app xx + ^pl- \j^^2 + \ in ^'( K )- (2-H) 

2.1. The case when a = 0. When <r = 0, ( |2.11| i becomes 

^ = ^-l^ + l-i^, iniy(R). (2.12) 

Since </? e and c - ^ we know that |c — is bounded away from 0. Hence 

from the standard local regularity theory to elliptic equation we see that <p £ C°° (M) and 
so is 77. Therefore in this case all solitary waves are smooth. 

As for the existence, we may multiply ( |2.12| i by <j) x and integrate on (—00, x] to get 

,„2 P 2 {c~p- A 1 ){c-p-A 2 ) 

ip x = r := &{<p), (2.13) 

c(c - ip) 

where 

A 1= ~ A + f^ , A 2 = - A ~f^ (2,4, 

are the two roots of the equation y 2 + Ay —1 = 0. Since A > 0, we know A\ > > A 2 . 
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From the decay property of <f> at infinity we know that a necessary condition for the 
existence is c > Ai or c < A 2 . But one may prove further that 

Theorem 2.5. When a = 0, ( |2.1| i admits a solitary wave solution if and only if 

c> A x or c<A 2 . (2.15) 

All solitary waves are smooth in this case. 

Proof. The regularity is discussed as above. So we will just focus on the existence part. 
If c = A\, then ( |2.13| l becomes 

,„2 -^{Ax- A 2 -ip) 

** = AMi-f) 1(V0 ' ( } 

Hence we see that <p(x) < near — oo. Because f(x) —> as x — > — oo, there is some 
Xq sufficiently large negative so that ip(xo) = — e < 0, with e sufficiently small, and 
<Px{xq) < 0. From standard ODE theory, we can generate a unique local solution (p(x) on 
[xq — L, xq + L] for some L > 0. Since A\ > > A 2 , we have 

p 2 [-3^ 2 + (6Ai - 2A 2 )p - 3Ai(Ai - A 2 )] 



V(Ai-A 2 -¥>) 



<0, 



(Ax - 

(2.17) 

for if < 0. Therefore Gi(ip) decreases for < 0. Because ip x (xo) < 0, if decreases 
near xq, so G\(<p) increases near xq. Hence from ( |2.16| l, <f x decreases near xq, and then 
ip and f x both decreases on [xq — L, xq + L]. Since JG\ {if) is locally Lipschitz in if for 
if < 0, we can easily continue the local solution to all of M and obtain that 0(x) — > — oo 
as x — ?> oo, which fails to be in H 1 (R). Thus there is no solitary wave in this case. 

Similarly we have that when c = A 2 there is no solitary wave. Therefore the theorem 
is proved. □ 

2.2. The case when a ^ 0. In this case we can rewrite ( |2.1 \\ as 

O--) 2 ) =fl- 2 -^^f+-f"-- + -^, inP'(K). (2.18) 

The following lemma deals with the regularity of the solitary waves. The idea is inspired 
by the study of the traveling waves of Camassa-Holm equation |33|. 

Lemma 2.6. Let a ^ and (if, rf) be a solitary wave of \2.\\ . Then 

(p- -) E C j (RV^Cc/o-)) , for k> 2 j . (2.19) 

Therefore 

ip£ C°° {®\<p-\c/a)) . (2.20) 

Proof. From Proposition 2.4 we know that c ^ and if ^ c and thus satisfies ( |2.18[ ). 
Let v = ip — — and denote 

= 3 /^\ 2 _ 2(c + A) /i, + e \ _ 1 
a \ cr J a \ a ) a 

So r(u) is a polynomial in v. From the fact that ip — c ^ we know that 

C — L c-i)^0. (2.21) 
fj 
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Then v satisfies 

-2 



(v ) X x = v x + r(v) H c - v 

cr V a 

From the assumption we know that (v 2 ) xx 6 Lj oc (R). Hence (v 2 ) x is absolutely continu- 
ous and hence 

v 2 G C^R), and then v G C 1 (R\u _1 (0)) . 
So from g2j) and that u + f G H 1 (M) C C(M) we know 

^-^c-vj GC , (R)nC 1 (M\u- 1 (0)). 
Moreover, 

(« fe ) ra = (^"V)* = ^ (v k -\v 2 ) x ) x 
= k(k-2)v k - 2 v 2 x + ^v k - 2 (v 2 ) xx 



k{k-2)v k - 2 v 2 x + k -v k - 2 



2 , , c 2 fa-1 

v x + rlv) H c — v 

a \ a 



kc 2 t _o ( a —\ 



= Hfc-- »"-X|-/"Mt)) + -»" c-« . (2.22) 

\ 2 J 2 2d y cr y 

For fc = 3, the right-hand side of ( |2.22[ ) is in L i 1 oc (]R). Thus we deduce that 

v 3 G C X (M). 

For fc > 4we see that ( |2.22| > implies 

(« fc ) M = 4 ( fc - 2) U ^ [(« 2 )J 2 +r fc " 2 H«)+^^- 2 (—c-,)" 2 G C(R). 

Therefore w fc € C 2 (R) for fc > 4. 

For fc > 8 we know from the above that 

v\ v k ~\ v k - 2 ,v k - 2 r{v) G C 2 (R), and v k - 2 (^^c -Jj 'eC 2 (RVr^O)) • 



Moreover we have 



v k - 2 v 2 x = \{v% ] ^ A {v k -%EC\M.). 



Hence from ( |2.22| > we conclude that 

v k G C 3 (MX-y-^O)) , fc>8. 

Applying the same argument to higher values of fc we prove that v k G C J (M\w _1 (0)) 
for fc > , and hence gig) . □ 

Denote 5 = min{a; : f(x) = c/a} (if tp ^ c/a for all x then let x — +00), then 
.t < +00. From Lemma 2.6 a solitary wave 99 is smooth on (—00, x) and hence \2.\ \\ 
holds pointwise on (—00, x). Therefore we may multiply by ip x and integrate on (—00, x] 
for x < x to get 

^= y2(c 7~y c "V A2) (2 - 23) 

(c - ¥j)(c-ct^) 
where Ai and A 2 are defined in ( |2.14| i. 



Applying the similar arguments as introduced in 11331 we make the following conclu- 
sions. 

1. When ip approaches a simple zero m = c — Ai or m = c — A2 of F(tp) so that 
F(m) = and F'(m) ^ 0. The the solution ip of ( |2.23| ) satisfies 

ip x = (ip — m)F' (m) + 0((ip — to) 2 ) as ip m, 

where / = 0(g) as x — >■ a means that \f(x)/g(x) \ is bounded in some interval [a — e, a+e] 
with e > 0. Hence 

ip(x) = m + -(x - x ) 2 F'(m) + 0((x - x ) 4 ) as x -> x , (2.24) 

where </?(x ) = to. 

2. If F(<p) has a double zero at ip = 0, so that F'(0) = 0, F"(0) > 0, then 

<f? x =<f?F"(p)+0(<p a ) as p-X). 

We get 

if(x) ~ aexp ^— ^^^"(O)^ as a; — > 00 (2.25) 

for some constant a. thus 93 — > exponentially as x — > 00. 

3. If ip approaches a simple pole ip(xo) — c/cr of i 7 '(</?) (when cr 7^ 1). Then 

- =/3|x-a; | 2/3 + 0((a;-Xo) 4/3 ) as x -J- x , (2.26) 




(2.27) 



for some constant f3. In particular, when F(tp) has a pole, the solution ip has a cusp. 

4. Peaked solitary waves occur when tp suddenly changes direction: ip x M> —<p x ac- 
cording to ( |2.23[ ). 

Now we give the following theorem on the existence of solitary waves of ( |2.1| i for 

a ^ 0. 

Theorem 2.7. For cr 7^ 0, we /zave 

(1) TfO < cr < 1, a solitary wave {ip, rj) of \2A) exists if and only if condition \2.\5\ 
holds. 

If c > A\ then ip > and max ie K f{x) — c — Ai. If c < Ai then ip < ant/ 
min xe K ^>(af) = c — A 2 . 

(2) //cr < 0, f/zen 

• if c> A\ then there is a smooth solitary wave ip > w/f/z max^gR = 
c — ^4i, ant/ an anticusped solitary wave (the solution profile has a cusp 
pointing downward) ip < with min^gR tp{x) = c/cr; 

• if c < A2 then there is a smooth solitary wave ip < with min xS R ip(x) — 
c — A2, and a cusped solitary wave ip > with max^gR ip(x) — c/a; 

• if c — A\ then there is an anticusped solitary wave ip < with min x gR ip(x) = 
c/a; 

• if c = A2 then there is a cusped solitary wave ip > with max^gR ip{x) = 
c/a. 

(3) If a > I, a solitary wave exists if and only if c satisfies ( |2.15| >. If c > A\ then 
ip > 0. If c < A 2 then ip < 0. Moreover, 
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• If A\ < c < — ztAi, then the solitary waves are smooth and unique up to 
translation with max ie i <p(x) = c — A\; 

• Ifc = — 3j Ai then the solitary wave is peaked with maXj.gR ip (x) = c—Ai = 
c/a; 

• Ifc > -^^A\ then the solitary waves are cusped with max l6 g f{x) = c/a; 

• // — £r-<4.2 < c < A 2 then the solitary waves are smooth and unique up to 
translation with min^R fix) = c — Ai; 

• If c= ~^z\A 2 then the solitary wave is antipeaked (the solution profile has a 
peak pointing downward) with min^R <p>(x) = c — A2 = c/a; 

• If c < 5~x^2 then the solitary waves are anticusped with min^gR ip(x) = 
c/a. 

Moreover, each kind of the above solitary waves is unique and even up to translations. 
When c > A\ or c < Ai, all solitary waves decay exponentially to zero at infinity. 

Proof. First from ( |2.23| l and the decay of <p(x) at infinity we know that a necessary condi- 
tion for the existence of solitary wave is that c > Ai or c < A 2 . 
If c — Ai then ( |2.23| > becomes 

2 -ip 3 (A 1 -A 2 -tp) . , . 
(^1 - WKM - vy) 



Hence we see that <p(x) < near —00. Similarly as in the proof of Theorem 2.5 



find some xq sufficiently large negative with <p(%o) = — e < and tp x (xo) < 0, and we 
can construct a unique local solution (p(x) on [xo — L-, x$ + L] for some L > 0. 

If a < 0, we see that A ^L av is decreasing when ip < 0. Together with ( |2.17| i we see that 
Fi((f) decreases for ip < 0. Because (p x (%o) < 0, <p decreases near xo, so Fi((p) increases 
near xq. Hence from ( |2.28| >, ip x decreases near xo, and then ip and <p x both decreases on 
[xq — L, xq + L], Since y/ F\(ip) is locally Lipschitz in tp for A\/a < ip < 0, we can 
easily continue the local solution to (—00, Xq — L] with ip(x) ~ > as x — > —00. As for 
x > xq + L, we can solve the initial valued problem 

ip(x + L) = p(x Q + L) 

all the way until ip = A\/a, which is a simple pole of Fi(ijj). From ( |2.26[ ) and ( |2.27| i we 
know that we can construct an anticusped solution with a cusp singularity at ip = A\/u = 
c/a. 

If a > 0, a direct computation shows that 

F[(<p) < 0, for <p < 0. 

Therefore the same argument indicates that <p{x) — > —00 as x — > +00, which fails to be 
in iJ 1 (E). Hence in this case there is no solitary wave. 

Similarly we conclude that when c = A 2 , there is no solitary wave when a > 0. When 
a < 0, there is a solitary wave with a cusp of height c/a. 

Now we consider c > A\ or c < A%. Again we will only look at c > Ay. The other 
case c < A 2 can be handled in a very similar way. From ( |2.23[ ) we see that p can not 
oscillate around zero near infinity. Let us consider the following two cases. 

Case 1. tp(x) > near —00. Then there is some x sufficiently large negative so that 
ip(xo) — e > 0, with e sufficiently small, and ip x (xo) > 0. 
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(i) When a < 1, \J F(<p) is locally Lipschitz in ip for < <p < c — Ai. Hence there is 
a local solution to 

{ <Px= V F (V)> 

on [xo — L,xq + L] for some L > 0. Therefore from ( |2.24| i and ( |2.25[ ) we see that in 
this case we can obtain a smooth solitary wave with maximum height ip — c — A\ and an 
exponential decay to zero at infinity 



ip(x) = O l^cxp |- Vc c C |s|JJ as |x|->oo. (2.29) 

(ii) When a > 1, yj F(ip) is locally Lipschitz in <p for < ip < c/a. Thus if c — A\ < 
c/a, i.e., c < -^tA\, it becomes the same as (i) and hence we obtain smooth solitary 
waves with exponential decay. 

If c — Ax = c/a then the smooth solution can be constructed until ip = c — A\ = c/a. 
However at <p = c — A\ = c/a it can make a sudden turn and so give rise to a peak. Since 
ip = is still a double zero of F(ip), we still have the exponential decay here. 

Lastly if c — Ai > c/a, then ip = c/a becomes a pole of F{tp). Hence from ( |2.26| l 
and ( |2.27| i we see that we obtain a solitary wave with a cusp at <p = c/a and decays 
exponentially. 

Case 2. <p{x) < near — oo. In this case we are solving 

f (p x = -y/F(<p), 
\ ip(x Q ) = -6 

for some xq sufficiently large negative and e > sufficiently small. 

When a > we see that F'(ip) < for ip < 0. Thus in this case there is no solitary 
wave. 

When a < 0, ip = c/a < is a pole of F(ip). Hence from similar argument as before, 
we obtain an anticusped solitary wave with minj.gR = c/a, which decays exponentially. 

Finally, from the standard ODE theory and the fact that the equation ( |2.11| > is invariant 
under the transformations x H> x + d for any constant d, and x H> —x, we conclude that 
the solitary waves obtained above are unique and even up to translations. 

□ 



Though there is no explicit expression for ip, and so rj in view of ( |2.10| i, as in 1581 . the 
effects of the traveling speed c on the function ip can be analyzed to provide some general 
description of its profile. Similarly to the case in [38 1 we have 



Proposition 2.8. Let c > A\ or c < A2, and ip is a smooth solitary wave of ( |2.1[ ) as 



obtained in Theorem 2.7 Then d c tp decays exponentially to zero at infinity and has at most 

2_ 
A' 



two zeros on K. In particular, if A\ < c < 4, the d c tp has exactly two zeros on . 



Proof. Again we only discuss the case c > A\. The other case c < A2 can be handled in 
the same way. 

Denote u> — d c ip. The exponential decay of ui can be inferred from ( |2.29| l. Since ip is 
unique and even up to translations, we may assume that p(0) = c — A\. Hence w(0) = 1 
and ui is even. Assume ui(xq) = for some xq > 0. Differentiating ( |2.23| l with respect to 
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c and evaluating at x = x we get 

2(p x u) x 



C — (TLp 
C — (Tip 



1 | {c-^) 2 +A{c- V )-l 
{c-tpf (c- <p){c- crip) 

1 



(c-ipy 



'r' 2 



> o, 



since c—aip > 0. Because (p x (xq) < 0, we see from the above inequality that lu x (xq) < 0. 
So lu is strictly decreasing near xq. It is then deduced from the continuity of oj that it has 
at most two zeros on R. 

If Ai < c < 4-, then from the decay estimate ( |2.29[ ) we see that cp decays faster at 
infinity as c gets larger, since 

/VSHM^rA 2- Ac >Q- 



2c Vc 2 + Ac - 1 



Hence < at infinity. Therefore oj has at least two zeros. Thus combining the above 
argument we proved that u>(x) has exactly two zeros ±xq in this case. □ 

Next we try to find an implicit formula for the peaked solitary waves. Let us consider 
only the case c > A%, Then from Theorem 2.7 we know that peaked solitary waves exist 
only when c = -^tA\. In this case we have 

2 <f 2 (c-A 2 -<p) 

fx = ■ 

c — ip 

Since <p is positive, even with respect to some xo and decreasing on (xq, oo), so for x > xo 
we have 



fx 




Hence from the separation of variables we get 



Let w = 1 



-{x - Xo) 



the above becomes 



dt 



-(x - x ) 



—A 2 

a 2 [cw — (c — A2)] (w — l)\/w 



dw 



A 2 



In 



cw — (c — A2) w 

CW - y/c- A 2 



dw 



+ s/c- A 2 



- In 


y/w — 1 




) 









Therefore we obtain an implicit formula for the peaked solitary waves. 

-\x - x \ = 



1 A 1 
1 A t 



A, 



In 



lew — v 1 c — A 2 



CW + y/c — A2 



- In 


y/w — 1 




1-- 


y/w + 1 


) 








w — 



A 2 

c-ip 



(2.30) 



Below is a figure of such a peaked solitary wave with xo = 0, A = 0, and c = a = 2. 
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3. Stability 

In this section, we want to discuss the stability of the smooth solitary waves of (2,1) . 
For fixed c and e > 0, we define the "e-tube" of a solitary wave cp c to be 

U t = {ueX : inf \\u- #.(• - s)\\ x < e}. (3.1) 

According to Theorem |2.7| the solitary waves for (2.1) travel with speeds proportional to 
their maximal heights. This consideration suggests that the appropriate notion of stability 
for the solitary waves is orbital stability: a wave starting close to a solitary wave should 
stay close, as long as it exists, to some translate of the solitary wave. The orbit of a solitary 
wave is the set of all its translates. 

Let us now discuss the appropriate notion of stability for the solitary waves of (2.1) . 

Definition 3.1. The solitary wave C of (2.1) is stable in X if for every e > 0, there exists 
a 5 > such that for any u £ Us, if u G C([0, T); X)for some < T < oo is a solution 
to (2.1) with u(0) = uq, then u(t) € U e for all t £ [0, T). Otherwise the solitary wave (p c 
is said to be unstable in X. 

As is discussed in (4), some solutions of (2.1) are defined globally in time (e.g. for 
< cr < 2 and inf^gR^o > —1, or the solitary waves constructed in Section |2| while 
other waves break in finite time. Note that by stability we mean that even if a solution 
which is initially close to a solitary wave blows up in a finite time, it will stay close to 
some translate of the solitary wave up to the breaking time. 

Our main theorem in this section is the following. 

Theorem 3.2. Let a < 1. All smooth solitary waves of (2.1) are stable. 



First from Theorem 2.7 we know that smooth solitary waves exist only when c > A\ or 
c < A 2 . For convenience we assume c > Ax- 

The special case a — 1 is settled in [38|. We will show that for general a < 1, the 
problem can be analyzed using the method provided by Grillakis, Shatah, and Strauss |26l . 

Let us now make some functional analysis setup. Recalling the functionals E and F are 
well-defined on X, we may compute their Frechet derivatives as follows 

E' u = -U xx +u f F' u = \ u 2 - - auu xx + r} + \ri 2 - Au 
E' v = r], \ F^ = U + UT). 
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Using these notation we see that a solitary wave ip c of ( |2.1[ ) satisfies 

cE'(<p c )-F'{(p c ) =0. (3.2) 

Denote 

L c = -dx ((c - crtp)d x ) - 3(p + atp xx + c + A. 
Then the linearized operator H c : X — > X* of C-E" — F' at can be computed as 



H c = cE"(0 c )-F"(0 c ) 
Using ( |2.10| i we have 



-(1+7?) C-ip 



( L * — \ 

C — (f 



(33) 



\ c-cp 

\ c — tp / 

We see easily that H c is self-adjoint and bounded from below, i.e., H c > al for some 
constant a and / is the identity operator. 

The next lemma states some spectral properties about H c . 

Lemma 3.3. Let c > A\ and (p c be a smooth solitary wave of ( |2.1| >. The spectrum of H c 
satisfies the following properties. 

(1) The essential spectrum of H c is positive and bounded away from zero. 

(2) The kernel of H c is spanned by d x (p c . 

(3) H c has exactly one negative simple eigenvalue Xi corresponding to eigenfunction 
X= (X,j"). 



Proof. The proof is inspired by Lemma 3.2 in |f38|. The details are as follows. 

(1) Since tp,ip x and <p xx all decay exponentially at infinity, it follows from Weyl's 
essential spectrum theorem that the essential spectrum of H c is the same as that of 
its asymptotic operator Ooo as |x| — > oo, where 



~cd xx +c + A -1 
-1 c 



When c > Ax, we have c 2 + Ac — 1 > 0. Hence there is some constant 5 
8(c, A) e (0, 1) such that 



2\iPu\ < 2(1 - 8)y/c{c + A)\ij)uj\ 
<(1-S) [{c + A)iP 2 + clo 2 ] 

for any $ = (ip,u) E H^R) x L 2 (R). Therefore 

(Ooo^, if T ) L 2 xL 2 = [ [cipl + (c + A)ij 2 - 2^w + cw 2 ] dx 

[6(c + A)iP 2 + Sclo 2 ] dx > dc\\ip\\ 2 L 2 xL 2. 



> 



Hence Ooo is positive when c > At, and then the essential spectrum of H c is 
[do, oo) for some do > and there are finitely many eigenvalues located to the left 
of a . 
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(2) Ifip — (ip, uj) is an eigenfunction of H c corresponding to the eigenvalue zero, then 



-d x ((c - a<p)ip x ) + (Sip + aip xx + c + A)ip u = 0, 

c — ip 

ip + (c — ip)ui = 0. 

c — p 



From the second equation we get ui = r^z^n ■ Hence the first equation can be 
expressed as a zero eigenvalue problem for JC C : H 1 — > H^ 1 : 



JC c \p := -d x ((c - aip)tp x ) + ( Sip + atp xx + c + A - ^_ ^ 3 ) 4> = 0. (3.4) 



We now use the fact that p(x), p x (x), ip xx (x) ~ * exponentially fast as |a;| — > 
oo while c—cr<p is positive and bounded away from zero when p is smooth. Similar 
to 1181 . it follows that the spectral equation lC c tp = can be transformed by the 
Liouville substitution 



. - d \ 9(z) = (c-*ip(x)) 1 / i tp(x), 
o \Jc-oip(y) 



into 



£ c 9(z) := [ -d 2 z + q c (z) +C + A--) 6(z) = 0, 



c 



where 

g c (z) = Sip (x) + — ip xx (x) - tt^ H ^7 j-^. 

4 (c~p{x)) i c 8[c—<Tip(x)) 

Since g c — > exponentially as |z| — ► oo, we deduce that £ c : — > H^ 1 (R) 

is self-adjoint with essential spectrum [c + A — -, oo). Because c > Ai, we know 
c + A — i > 0. We may have finitely many eigenvalues of C c located to the 
left of c + A — -. The nth eigenvalue (in increasing order) has, up to a constant 
multiple, a unique eigenfunction with precisely (n — 1) zeros (see for example, 
[24 1 for details). Thus the operator JC C has the same spectral properties. 

Note that r/ x = (°-(p) 2 anc ' *H3 i m ply that fC c (ip x ) — 0. Since ip x has exactly 
one zero. Therefore the zero eigenvalue of JC C is simple, and there is exactly one 
negative eigenvalue while the rest of the spectrum is positive and bounded away 
from zero. Hence the zero eigenvalue of H c is simple and the kernel is spanned 
by <p x . 

(3) The operator H c is related to a quadratic form Q c (tp) with ip — (tp,oj) € X, 
which is defined as the coefficient of e 2 in the Taylor's expansion of cE((p c + 
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etp) — F(if c + etp) and is given by 
1 



=W0 



(c - <J^)i\) x + (—3(^5 + cr^Kx + c + A)^ 



2c 



+ (C - </?)t</ 



(c - cry>)V^ 



5- V) 



r 0- W 



c + A- 



dx+ 



dx 



Note that the quadratic form (ip) is related to the operator K c and G(if) is 
nonnegative. 

Let / be a nontrivial eigenfunction corresponding to the unique negative eigen- 
value of JC r - Then 



f)=Q { c 1] (f)<o 



Qdf, 



So H c has a negative eigenvalue, say, Ai < 0. Applying the min-max characteri- 
zation of eigenvalues to H c yields 

A2 = max 



*ex <2eX\{0} \H\x 
(£&J) = 



A 2 > 



Choosing -0 = (/, 0) leads to 

<2c(w) 



Qi 1) (g) + G(g,M ^ n 

(g,h)eX\{0} \\{9,h)\\ x 
((f,0),(g,h))=0 



mm 

a e x\{o} 

((/,0),£3) = 



A" 



The last inequality is due to Qc (<?) > for all g such that (g, = and that 
G(g, h) > 0. Therefore Ai is simple. Denote the corresponding eigenfunction 
by x — (XiA*)- From the result in (2) we see that A2 = is also simple. This 
completes the proof of the lemma. 

□ 

Remark 3.4. (i) Notice that the above lemma applies to all smooth solitary waves of j2.1\ 
without the restriction that a < 1. 

Under the assumption c < A2, we can consider the operator 

c \ / c \ 



( L c 



H r = - 



\ 



c — if 
C — if 



( -L c 



V 



c — if 



c — ip 

— C + if 



c — ip 

In this setting we have for smooth solitary wave that c < c — A2 < < and c — aip 
is bounded away from zero. By a similar argument, all properties of H c in Lemma 3.3 are 
still valid. 
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(ii) We will apply the method of Grillakis-Shatah-Strauss [26] to establish the stability 
of smooth solitary waves. However our problem does not exactly fall into the frame work 
there since the operator J is not onto. But in fact the invertibility of J is only needed to get 
instability and is not required for stability (see Sections 3 and 4 in [26 1 for more details). 
Hence the argument in Section 3 of [26 1 can still be used here. 



Let (p c = (ip, -q) be a solitary wave of \2.\\ . Consider the following scalar function 



d(c) 



cE(0 c ) - F(0 C ) if c> A u 
F(0 C ) - cE(<p c ) if c<A 2 . 



(3.5) 



The next lemma shows that for a < 1, , c > Ay or c < A 2 and ip c = (ip, rj) G X being a 
smooth solitary wave of \2.\\ , d(c) is convex in c. 

Lemma 3.5. Assume a < 1, c > Ay or c < A% and tp c = (ip, rj) is a smooth solitary wave 
of§n\. Then d"(c) > 0. 

Proof. Consider first c > A\. Differentiating d(c) with respect to c and then applying ( |3.2| i 
we obtain 



d'(c) = (cE'(0 c ) - F'(p c ),d c <p c ) + E(p c ) - E{<p c ). 
In view of the even symmetry of <p, it follows from ( |2.10[ ) and ( |2.23| ) that 

■2 



(3.6) 



d'(c) = E{0 C ) 



<p£ + <r + 



dx 



(c-tp) 2 
{c - tp - Ay)(c - tp - A 2 ) 

(C- <P)(C- (Tip) 



+ 1 



{c-pf 



dx. 



Recall that < <p{x) < c — A\ and (p'(x) < on (0, oo) when c > A\. We have from 
( |2T23] > that 



d'(c) 



(C- tp)(c~ (Tip) 



(c-cp-A^c-tp-Az) 
(c-tp-A 1 )(c-tp-A 2 ) , , , 1 



(c - tp)(c - atp) (c-ip) 2 
Introducing a change of variable y — c — ip(x) the above becomes 

"(y- A 1 )(y-A 2 ) , 



dx. 



(3.7) 



d'(c) 



[c~y)\ 



y [(1 - cr)c + cry] 
(y-^i)(y-A 2 ) 



2/ [(1 - cr)c + ay] 



dx. 



Differentiating the above with respect to c we have 



d"(c)=0 + dMc-y) s 



+ del {c-y)\ 



h(y) dy- 



\y-A 1 )(y-A 2 ) 
y [(1 - cr)c + ay] 

I y [(1 - (r)c + oy\ 
(y-A 1 )(y-A 2 ) 

h{y) dy. 



dy 



dy 



Since ip is smooth, we have that c — ap > 0. Hence (1 — a)c + cry = c — cr<p > 0. 
Moreover, we know that A 2 < < A\ < y < c. Further explicit computation shows that 



T( s = (V- M){y~ A 2 ) f [(l-a)c + ay}+y 
1{V) V v[Q--<r)c + w] ' V 2[(l-a)c + ay] 
Here we don't need to assume a < 1. If a < 1, we have 

(y 2 + l) [(l-a)c + ay + ±(l-a)(c-y)] 

h{y) - 



> 0. 



> 0. 



□ 



y V (y- A i)(y- ^M(i - a)c + o-y 

Therefore cf'(c) > 0. 

The other case c < A 2 can be handled in a very similar way and hence we omit it 

The next lemma can be obtained by exactly the same proof as in ll26l . 

Lemma 3.6. Let a < 1 and <p c = (ip, if) be a solitary wave of (2,1) . There exist e > 
and a unique C l map s : U e — > K 5mc/i that for every u £ U e and r£l 

(1) 

+ s(u)),ip x ) = 0, 

(2) 

s(u(- + r)) = s(u) — r. 



By the spectrum analysis in Lemma 3.3 and the convexity of d(c) in Lemma 3.5 



can follow exactly the same idea as in [26 1 Theorem 3.3 to get 



Lemma 3.7. Let the assumptions of Lemma 3.5 hold. There exists a constant k — k(c) > 
such that 

(H c $),$ >k\\$\\ 2 x , (3.8) 
for all ip £ X satisfying (ip c , ip) — (cp' cl ip) = 0. 



The following lemma can be obtained directly from Lemma [3T6| and Lemma |377| 
Lemma 3.8. Let the assumptions of Lemma \3.5\ hold. There exists an e > such that 

F{0 c )-F{u)>^\\u{- + s{u))-(p c \\ 2 x , 
for u £ U e satisfying E(u) = E(tp c ). 



Proof of Theorem \3. 2\ In view of Lemma 3.3 and Lemma 3.8 the result of theorem is then 
a direct consequence of Theorem 3.5 in [26 1. □ 



4. Global solutions when a = 

In Ej, the authors established a blow-up criterion for a ^ (cf. Theorem 3.3 in (4|). In 
fact, the restriction of a ^ can be removed using the same argument and hence we get 

Theorem 4.1. Let (u, p) be the solution of \l.2\ with initial data (u$, po — 1) £ H S (M.) x 
H S ~ 1 (M), s > 3/2, andT the maximal time of existence. Then 

T<^ [ K(r)||ii-dr = oo. (4.1) 



The wave-breaking phenomena for system ( |1.2) when a ^ was discussed in details 
in |4) . Here we show that when a = the solutions constructed in Theorem |2.2| are 
global-in-time. 
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Theorem 4.2. Let a = 0. If(uo,po — 1) £ H s x H s \ s > 3/2, f/zen f/zere exz'sfs a 
unique solution {u,p-l) of in C([0, oo);H s x H"- 1 ) n C^flO, oo); i? 8 " 1 x # s ~ 2 ) 
wzf/z (u, p) |t=o = ( u o, Po)- Moreover, the solution depends continuously on the initial data 
and the Hamiltonian H\ is independent of the existence time. 

As discussed in [4j, system \l.2\ has two associated characteristics q and q given by the 
following initial-value problems 

{ 0<i<T, (42) 

[ g(0,x) =x, x e M, 

— =au{t,q), 0<t<T, (43) 

g(0,a;)=a; ) a: G E, 

where it € C 1 ([0,T),i?' a-1 ) is the first component of the solution (u, p) to \\.2\ with 
initial data (uo,Po) & H s x with s > 3/2 and T > is the maximal time of 

existence. When a = 0, the second one q becomes stationary. Thus we will perform the 
estimates along the first characteristics q. 

A direct calculation shows that for t > 0, x £ R 

q x {t, X ) = e Jtv x (T,q(T,x))dT > Q 

Hence q(t, •) : E — > M is a diffeomorphism of the line for each t £ [0, T). Hence the L°° 

norm of any function ■) g g [0,T) is preserved under q(t, ■) witht G [0,T), 
i.e., 

IKt,-)IU-(R) = |Kt,?(t,-))lli-(R), *e[0,T). (4.4) 

Similarly we have 

inf u(t,ar) = inf u(*,g(t,a;)), te[0,T), (4.5) 

supt;(i,x) = supt>(t,g(t,x)), ie[0,T). (4.6) 

When cr = 0, we can rewrite system jl.2\ as 

u t + d x p * (-Au + |u 2 + \p 2 ) = 0, 
Pt + (Hx = 0, 

where p(x) is defined in ( |2.5| > 

The following lemma is needed in candying out the estimates along the "extremal" char- 
acteristics. 

Lemma 4.3. (M)LetT > Oandv £ C 1 ([0, T); H 2 (R)). Then for every t £ [0,T) there 
exists at least one point ((t) 6 I with 

m(t) := inf [v x {t,x)\ = ^ (*,£(*)) ■ 
The function m(t) is absolutely continuous on (0, T) with 

=«**(*,£(<)) o«(0,T). 

To prove Theorem |4.2| of global well-posedness of solutions, we need the following 
estimates for u x . 
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Lemma 4.4. Let a — and (u, p) be the solution of ( |4.7| > with initial data (v,q, po — 1) € 

_ff s (R) x JP _1 (K), s > 3/2, and T the maximal time of existence. Then 

sup u x (t,x) < svpu 0lX (x) + - ( suppl(x) + Cf)t, (4.8) 



1 



inf u x (t,x) > inf u 0tX (x) + - inf ^(x) - C£ i, (4.9) 



where the constants above are defined as follows. 



Ci = V^— IIKpo-i)!!*^, (4.io) 



C 2 = y]2 + Ct (4.11) 

Proof. The local well-posedness theorem and a density argument implies that it suffices 
to prove the desired estimates for s > 3. Thus we take s — 3 in the proof. Also we may 
assume that 

u # 0. (4.12) 
Otherwise the results become trivial. Since now s > 3, we have u <E Cq(M). Therefore 

inf ti x (t,x) < 0, supu x (t,x) > 0, ie[0,T). (4.13) 

Differentiating the first equation of ( |4.7[ ) with respect to x and using the identity —d 2 p* 
f = f — p * / we obtain 

u tx = ^p 2 + '^u 2 + Adlp*u~p* (\u 2 + |p 2 J . (4.14) 



Using Lemma 4.3 and the fact that 

sup (t, x)] = - inf [-v x (t, x)} , 

we can consider m(i) and £(t) as follows 

m(t) :=«*(*,£(*)) = sup («,(*,«)), te[0,T). (4.15) 

x£R 

Hence 

u ra (t,e»)-0, a.e. ie[0,T). (4.16) 
Take the trajectory q(t,x) defined in ( |4.2[ ). Then we know that g(t, •) : M. —> R is a 
diffeomorphism for every i £ [0, T). Therefore there exists Xi(t) G R such that 

q(t, Xl (t))=at) te[0,T). (4.17) 

Now let 

C(*) = p(t,g(t,a! 1 )), te[o,r). (4.18) 

Therefore along this trajectory g(t, Xi) equation ( |4. 14| i and the second equation of ( |4.7| i 
become 

m'^^C 2 + /(*,<?(*, zi)), 

C'(f) = -Cm, (4.19) 
for i € [0, T), where ' denotes the derivative with respect to £ and /(t, g(t, x)) is given by 

/= ^u 2 + Ad 2 xP *u-p* (1^ + \p 2 )- (4-20) 



We first derive the upper and lower bounds for / for later use in getting the wave- 
breaking result. Using that d 2 p * u — p x * u x , we have 

3 3 1 1 

/ = + A Px * u x - ^P*u 2 - -p*l - p* {p - 1) - -p* (p- l) 2 

< ^u 2 +A\ Px *u x \-- + \p*(p-l)\. 

Since 

A\p x *u x \ < A\\p x \\ L 2\\u x \\ L 2 = -A\\u x \\ L 2 < - + ^A 2 \\u x \\ 2 L2 , (4.21) 

|p*(p-l)| < ||p|M|p-l|U' = <\ + \\\p-M\h, (4-22) 

u 2 < - [ (u 2 + u 2 x ) dx, (4.23) 
we obtain the upper bound of / 

f< \\\P- M\h + l\\u\\h + ^^\\U X \\ L 2 

< ^-||(«o,P0 - i)\\mxv = \cl (4.24) 

Now we turn to the lower bound of /. Similar as before, we get 

3 1 1 

-/ < A\p x * u x \ + -p * u 2 + - + \p * (p - 1)| + -p * {p - l) 2 . 



- \ + ^\WAh + \\W\\h + \\\p - i h 



< i + ^±^!|| (M05po _ l)||^ lxi2 = hoi (4-25) 



where we have used the inequality 



P*g 2 < \\\9 2 \\li = \h\\h- 

Combining ( |4.24| > and ( |4.25| l we obtain 

l/l < 1 + ^-||(«o, A) - 1)II^xl- (4-26) 

From ( |4.13| ) we know fh(t) > for t € [0, T). From the second equation of ( |4.19| l we 
obtain that 

((i) = ((0) e - rn(T)dT . (4.27) 

Hence 

| P (t j9 (t ia:i ))| = |C(t)|<|C(0)|. 

Therefore we have 

fh'(t) = k 2 (t) + / < k" 2 (0) + \C 2 < \ (suvp 2 {x) + C\ 



Integrating the above from over [0, t] we prove ( |4.8) , 

To obtain a lower bound for inf^gR u x (t 1 x), we use the similar idea. Consider the 
functions m(t) and as in Lemma [43] 

m{t) := u x (t, £(t)) = inf ( Ujs (t, a;)) , t G [0, T). (4.28) 
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Hence 

Ux» (*,€(*)) = a.e. te[0,T). (4.29) 
Again take the characteristics g(t, x) defined in j4-.2\ and choose x 2 (t) € K such that 

g(t,a; 2 (t))=£(i) te[0,T). (4.30) 

Let 

C(t)=p(t,q(t,X2)), te[0,T). (4.31) 
Hence along this trajectory q(t,X2) equation ( |4.14[ i and the second equation of j4.7\ be- 
come 

m\t)= l -C 2 + f(t,q{t,x 2 )), 

C(t) = -Cm. (4.32) 
Since m(t) > 0, we have from the second equation of the above that 

| P (t,«(t )a:a ))| = |C(t)|>|C(o)|. 

Then 

Integrating the above from over [0, t] we obtain ( |4.9| ). This completes the proof of Lemma 
Lemmal4~4l □ 



i'(t) > -C 2 (0) - -C\ > - iniplix) - C\ 



Proof of Theorem \4. 2| Combining Lemma 4.4 and Theorem 4. 1 we easily see that the local 



solution obtained in Theorem 2.2 can be extended to all of the interval [0, oo). □ 



References 

[1] R. Beals, D. Sattinger and J. Szmigielski, Multi-peakons and a theorem of Stieltjes, Inverse Problems 15 

(1999) , L1-L4. 

[2] T.B. Benjamin, J.L. Bona and J. J. Mahony, Model equations for long waves in nonlinear dispersive systems, 

Philos. Trans. Roy. Soc. London Ser. A 272 (1972), 47-78. 
[3] R. Camassa and D. Holm, An integrable shallow water equation with peaked solitons, Phys. Rev. Lett. 71 

(1993), 1661-1664. 

[4] R.M. Chen and Y. Liu, Wave-breaking and global existence for a generalized two-component Camassa- 

Holm system, Int. Math. Res. Not., (2010), in press. 
[5] A. Constantin, On the Cauchy problem for the periodic Camassa-Holm equation, J. Differential Equations, 

141(1997), 218-235. 

[6] A. Constantin, Global existence of solutions and breaking waves for a shallow water equation: A geometric 

approach, Ann. Inst. Fourier (Grenoble), 50 (2000), 321-362. 
[7] A. Constantin, On the blow-up of solutions of a periodic shallow water equation, J. Nonlinear Sci., 10 

(2000) , 391-399. 

[8] A. Constantin, The trajectories of particles in Stokes waves, Invent. Math., 166 (2006), 523-535. 
[9] A. Constantin and J. Escher, Wave breaking for nonlinear nonlocal shallow water equations, Acta Math., 
181(1998), 229-243. 

[10] A. Constantin and J. Escher, Well-posedness, global existence and blow-up phenomena for a periodic quasi- 
linear hyperbolic equation, Comm. Pure Appl. Math., 51 (1998), 475-504. 

[11] A. Constantin and R. Ivanov, On the integrable two-component Camassa-Holm shallow water system, Phys. 
Lett. A, 372 (2008), 7129-7132. 

[12] A. Constantin and D. Lannes, The hydrodynamical relevance of the Camassa-Holm and Degasperis-Procesi 
equations, Arch. Ration. Mech. Anal., 192 (2009), 165-186. 

[13] A. Constantin and H. P. McKean, A shallow water equation on the circle, Comm. Pure Appl. Math., 52 
(1999), 949-982. 

21 



[14] A. Constantin and R.I. Ivanov, On an integrable two-component Camassa-Holm shallow water system, 

Phys. Lett. A, 372 (2008), 7129-7132. 
[15] A. Constantin and L. Molinet, Orbital stability of solitary waves for a shallow water equation, Phys. D 157 

(2001) 75-89. 

[16] A. Constantin and W.A. Strauss, Stability ofpeakons, Comm. Pure Appl. Math., 53 (2000), 603-610. 
[17] A. Constantin and W.A. Strauss, Stability of the Camassa-Holm solitons, J. Nonlinear Sci. 12 (2002), 415— 
422. 

[18] A. Constantin and W.A. Strauss, Stability of a class of solitary waves in compressible elastic rods, Phys. 

Lett. A 270 (2000), 140-148. 
[19] H.-H. Dai, Model equations for nonlinear dispersive waves in a compressible Mooney-Rivlin rod, Acta 

Mech., 127 (1998), 193-207. 
[20] PG. Drazin, Solitons, London Mathematical Society Lecture Note Series, Cambridge University Press, 

Cambridge, 1983. 

[21] H. Dullin, G. Gottwald and D. Holm, On asymptotically equivalent shallow water wave equations, Physica 
D, 190 (2004) 1-14. 

[22] H. Dullin, G. Gottwald and D. Holm, Camassa-Holm, Korteweg-de Vries-5 and other asymptotically equiv- 
alent equations for shallow water waves, Fluid Dyn. Res., 33 (2003) 73-95. 

[23] H. Dullin, G. Gottwald and D. Holm, An integrable shallow water equation with linear and nonlinear 
dispersion, Phys. Rev. Lett., 87, (2001) 194501-04. 

[24] N. Dunford and J.T. Schwartz, Linear Operators. Part II: Spectral Theory, Selfadjoint Operators in Hilbert 
Spaces, Wiley, New York, 1988. 

[25] A. Fokas and B. Fuchssteiner, Symplectic structures, their B'dcklund transformation and hereditary symme- 
tries, Phys. D, 4 (1981), 47-66. 

[26] M. Grillakis, J. Shatah, and W. Strauss, Stability theory of solitary waves in the presence of symmetry, J. 
Funct. Anal. 74 (1987), 160-197. 

[27] G. Gui and Y. Liu, On the Cauchy problem for the two-component Camassa-Holm system, Math. Z., (2010), 
DOI: 10.1007/s00209-009-0660-2. 

[28] D. Holm, L Naraigh and C. Tronci, Singular solutions of a modified two-component Camassa-Holm equa- 
ripn larXiv:0809.2538l 

[29] R. Ivanov, Two-component integrable systems modelling shallow water waves.The constant vorticity case, 

Wave Motion 46 (2009), 389-396. 
[30] R.S. Johnson, Camassa-Holm, Korteweg-de Vries and related models for water waves, J. Fluid Mech. 455 

(2002) , 63-82. 

[31] R.S. Johnson, The Camassa-Holm equation for water waves moving over a shear flow. Fluid Dynam. Res. 
33 (2003), 97-111. 

[32] C. Kenig, G. Ponce and L. Vega, Well-posedness and scattering results for the generalized Korteweg-de 
Vries equation via the contraction principle, Comm. Pure Appl. Math., 46 (1993), 527-620. 

[33] J. Lenells, Traveling wave solutions of the Camassa-Holm equation, J. Differential Equations, 217 (2005), 
393^30. 

[34] J. Lenells, Traveling waves in compressible elastic rods, Discrete Contin. Dyn. Syst., 6 (2006), 151-167. 
[35] J. Lenells, Stability of periodic peakons, Int. Math. Res. Not., 10 (2004), 485^199. 

[36] J. Lenells, A variational approach to the stability of periodic peakons, J. Nonlinear Math Phys., 11 (2) 
(2004) 151-163. 

[37] Y. Li and P. Olver, Convergence of solitary-wave solutions in a perturbed bi-Hamiltonian dynamical system. 

I. Compactons and peakons Discrete Contin. Dyn. Syst., 3 (1997), 419^132. 
[38] Y. Liu and P. Zhang, Stability of solitary waves and wave-breaking phenomena for the two-component 

Camassa-Holm system, Int. Math. Res. Not., 211 (2009), 41 pages. 
[39] H. P. Mckean, Integrable systems and algebraic curves, Global Analysis, Springer Lecture Notes in Mathe- 
matics, Vol. 755, Springer, Berlin, 1979, 83-200. 867-874. 
[40] O. Mustafa, On smooth traveling waves of an integrable two-component Camassa-Holm shallow water 

system, Wave Motion 46 (2009), 397^102. 
[41] G.A. Nariboli, Nonlinear longitudinal dispersive waves in elastic rods, J. Math. Phys. Sci., 4 (1970), 64-73. 
[42] A. Shabat, L. Martinez Alonso, On the prolongation of a hierarchy of hydrodynamic chains, Proceedings 

of the NATO advanced research workshop, Cadiz, Spain 2002, NATO Science Series, Kluwer Academic 

Publishers, Dordrecht, (2004), 263-280. 
[43] W.A. Strauss, Nonlinear Wave Equtions, Conf. Board Math. Sci., Vol 73, American Mathematical Society, 

Providence, RI, 1989. 
[44] J. F. Toland, Stokes waves, Topol. Methods Nonlinear Anal., 7 (1996), 1—48. 



22 



[45] G.B. Whitham, Linear and Nonlinear Waves, John Wiley & Sons, New York, 1980. 
Robin Ming Chen 

School of Mathematics, University of Minnesota, Minneapolis, MN 55455 
E-mail address: chenm@math . umn . edu 

Yue Liu 

University of Texas at Arlington, Department of Mathematics, Arlington, TX 76019-0408 
E-mail address: yliu@uta.edu 

Zhijun Qiao 

Department of Mathematics, University of Texas-Pan American, Edinburg, TX 78539 
E-mail address: qiao@utpa.edu 



23 



